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We investigate the spin tunnelling of spin-1 Bose-Einstein condensates in a linearly sweeping 
magnetic field with a mean-field treatment. We focus on the two typical alkali Bose atoms Rh 
and '^^Na condensates and study their tunnelling dynamics according to different sweeping rates of 
external magnetic fields. In the adiabatic (i.e., slowly sweeping) and sudden (i.e., fast sweeping) 
limits, no tunnelling is observed. For the case of moderate sweeping rates, the tunnelling dynamics is 
found to be very sensitive on the sweeping rates with showing a chaotic-like tunnelling regime. With 
magnifying the regime, however, we find interestedly that the plottings become resolvable under a 
resolution of 10"'* G/s where the tunnelling probability with respect to the sweeping rate shows 
a regular periodic-like pattern. Moreover, a conserved quantity standing for the magnetization in 
experiments is found can dramatically affect the above picture of the spin tunnelling. Theoretically 
we have given a reasonable interpretation to the above findings and hope our studies would bring 
more attention to spin tunnelling experimentally. 

PACS numbers: 03.75.Mn, 03.75.Kk 
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I. INTRODUCTION 

Bose-Einstein condensation (BEG) has been one of the 
most active topics in physics for over a decade, and yet 
interest in this field remains impressively high. One of 
the hallmarks of BEG in dilute atomic gases is the rel- 
atively weak and well-characterized interatomic interac- 
tions. The vast majority of theoretical and experimen- 
tal work has involved single component systems, using 
magnetic traps confining just one Zeeman sublevel in the 
ground state hyperfine manifold, including the BEG-BGS 
crossover [l|, , quantized vortices 0, 0) Si j condensates 
in optical lattices @ , and low-dimensional quantum gases 
0, [M] ■ An impotent frontier in BEG research is the exten- 
sion to multicomponent systems, which provides a unique 
opportunity for exploring coupled, interacting quantum 
fluids. In particular, atomic BEGs with internal quantum 
structures, some experiments have observed spin proper- 
ties of F = 1 and F = 2 condensates i, [Ulli!, Ill, til , 
using a far-off resonant optical trap to liberate the in- 
ternal spin degrees of freedom. Even F = 3 bosons are 
also investigated in a present theoretical work For 
atoms in the F = 1 ground state manifold, the presence 
of Zeeman deg enerac y and spin-dependent atom-atom in- 
teractions [l3,[l^[la[l3j[I3j[H leads to interesting con- 
densate spin dynamics, especially spin-1 system with its 
relatively simple internal structure. Many literatures has 
been devoted to spin mixing and spin domain which had 
been observed in experiments. 

In this article, we investigate spin tunnelling for spin- 
1 BEG with a mean-field description. Unlike all pre- 
vious studies of the fixed external magnetic fields both 
in theory and experiment (e.g., see Refs [H, [H, H^j), 
we highlight the important role of an external magnetic 
field that is now set to be linearly varying with time. 
We focus on the two typical alkali Bose atoms and 



'^^Na condensates and study their tunnelling dynamics 
according to different sweeping rates of external mag- 
netic fields. We also pay much attention to a conserved 
quantity, m, standing for magnetization, and find that 
this quantity can dramatically affects tunnelling dynam- 
ics for both Rb and '^^Na atom system. 

Our paper is organized as follows. Sec. II introduces our 
model. In Sec. Ill we demonstrate our numerical simula- 
tions on the irregular spin tunnelling of alkali Bose atoms 
^'^Rb and "^^Na condensates respectively. In Sec. VI, we 
present a theoretical interpretation to the above findings 
with the help of both analytical deductions and phase 
space analysis. Sec.V is our conclusion. 



II. THE MODEL 

In an external magnetic field, spin-1 Bose-Einstein con- 
densate (BEG) is described with the following Hamilto- 
nian [l^ 

H = J dr[V'l(-|^V2+«+F0^.+|vlV'jV',V'.+yV'iV'l(i^7)y(^^7 

(1) 

where repeated indices are summed and ipjir) (^pi) is the 
field operator that creates (annihilates) an atom in the 
i-th hyperfine state {\F = 1, i = -|-1, 0, — 1), hereafter \i)) 
at location r. M is the mass of an atom. Interaction 
terms with coefficients cq and C2 describe, respectively, 
elastic collisions of spin-1 atoms, expressed in terms of 
the scattering length ag (a2) for two spin-1 atoms in the 
combined symmetric channel of total spin (2), cq — 
AirTi^ (ao + 202) /3M and C2 = 47rS.^ (as - oq) /3M. qq is 
not spin concerned. 02 is spin concerned. F^=x,y,z are 
spin-1 matrices. Assuming the external magnetic field B 
to be along the quantization axis (z), the Zeeman shift 
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on an atom in state \i) becomes 
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where Eh j is the hyperfine spUtting and gi is the Lande 
g factor for an atom with nuclear spin /. fij is the nu- 
clear magneton and ^ = (gi^ijB + gj^sB) /Ehf with gj 
representing Lande g factor for a valence electron with a 
total angular momentum J. /is is the Bohr magneton. 

At near-zero temperature and when the total number 
of condensed atoms (N) is large, the system can be well 
described in the mean-field approximation. For isotropic 
Bose gas, under the mean-field method and single model 
approximation, the operators can be substituted with c 
numbers ipi = ai4>{r) where ai correspond to the prob- 
ability amplitudes of atoms on i-th hyperfine state. By 
setting ai = ^fsle^^^ , the system can be described by the 
following classical Hamiltionian system [2^, 
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S+si+£;oSo+-B-S-i-c[(l-si-s_i)^-|-4sis_i-4(l- 

(4) 

where Q = Q\ + Q-\ — 29o and c = C2 J dr |(/)(r)| . Using 
canonically conjugate transformation, Hmf can be trans- 
fered into the following compact classical Hamitonian (up 
to a trivial constant) 



= 2cso[(l - so) + V(l-5o)2-m2 cos0] + 5{1 - sq), 

(5) 

and equations of motions for canonically conjugate vari- 
ables So, d are 



So 



4c r— 



so)2 — m? sin( 



2(5 



4c , , 4c (l~so)(l-2so)~m2 
-(l-2so) + ^ 



(6) 



cos 6*, 
(7) 

where m = s_i — siis conserved and denoted as magne- 
tization, and 5 — {E+ + E- — 2Eq)/2. Fig.l shows the 
relationship between 5 and the external magnetic field B. 



III. IRREGULAR SPIN TUNNELLING 

As one of spin dynamics problems, spin tunnelling is 
always interesting to theoretical and laboratorial inves- 
tigations. In this section, we study spin tunnelling for 
the spin-1 BEC systems in a sweeping magnetic field. 
In our study, the magnetic field varies in time linearly 
, i.e., ~ ai, from B —> —oo to B — > -|-oo. The sweep- 
ing is far away from Feshbach resonance and ensure that 
the atom-atom interaction is almost not variety during 
the sweeping process. Numerically, the interval of the 




FIG. 1: The transformation of S, expressed in term of the 
.^e£iyamsljj^t^DTjia^g(|om in state |i) (i — -1-1,0,-1) in the 
external magnetic field B. (a) for Rb atom condensate, (b) 
for ^^Na atom condensate. 



magnetic fields is set as [—Bo,Bo] and the Bq is cho- 
sen larger enough so that the coupling between different 
components are safely ignored at beginning and ending. 
In this situation, the spin tunnelling probability can be 
well defined. We want to know the final value of sq ( i.e. 
) at _B — > +oo{i.e., Bq), suppose initially we have 
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at i? — > —oo{i.e., —Bq). We exploit Runge-Kutta 
^th _ t^th algorithm to numerically solve the coupled or- 
dinary differential equations ([6]) ([7]) for the parameters 
corresponding to Bose atoms Rh and '^^Na respectively. 



A. Atom Condensate 

Because ^"^ Rb atom condensate can be readily prepared 
in experiments and has been involved in many investiga- 
tions, we firstly discus it. For convenience and not losing 
generality, we set the initial probability of so as 0.5 for an 
example. Fig. 2 plots the final value of so at B ^ -l-oo, 
i.e. s^™°' for different sweeping rates. The above plot- 
ting suggests that our discussions on the spin tunnelling 
can be divided into three parts according to the values 
of the sweeping rates a. 

a) When a ^ 0, s^*"°' almost equal to initial con- 
dition So = 0.5, as if the system has not been changed. 
The plotting of tunnelling probability vs. sweeping rates 
is almost a line, which indicates that no tunnelling occurs 
after the external magnetic fields sweep slowly from neg- 
ative infinity to positive infinity. In this case, the system 
is believed to adiabatically change with the slowly sweep- 
ing magnetic field. When the magnetic field changes from 
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FIG. 2: The tunnelling probability of ^"^Rb for m = 
0.0, 0.1, 0.3 at c = —3.13 x 10"^'' J in the sweeping magnetic 
field. 



— oo to +CXD, S, the Zeeman shift, changes from an initial 
value to zero, then returns back to its initial value be- 
cause of its symmetric dependence on the field (shown 
in fig.l). Hence, the classical Hamiltonian ([5]) comes 
back to origin though the field has been changed dra- 
matically. Numerically solving Eqs.([6l [7]) (also see phase 
space fig. 9) we see that for a fixed magnetic field, the 
motions are periodic. Hence, in light of the adiabatic 
theory 21], if the sweeping rate of the external magnetic 
field is small compared to the frequencies of the instanta- 
neous periodic orbits, the the system undergoes adiabatic 



evolution. Therefore, Sq « sq. 

b) When a — > oo, s^™°' also tends to 0.5 (its inital 
value). The final value of sq oscillates around 0.5 and 
tends to a line. In this case, the sweeping rate is so quick, 
and the system ([5]) restore quickly. If the time of change 
the magnetic field is much shorter than the peroid of the 
motion of system. It is expected that there is no time for 
the system to give some response to the change of field. 
So no tunnelling phenomenon for very fast sweeping rate 
can also be well comprehended. 

c) The interesting phenomena emerge when a is mod- 
erate. For this case, we find s^™"' changes dramatically 
with respect to sweeping rates, which indicates spin tun- 
nelling occurs and the tunnelling probability is seemingly 
chaotic [23|. The spin tunnelling process can be shown by 
drawing the evolution of sq with respect to instantaneous 
magnetic fields B. In fig. 3, we plot the temporal evolution 
of So for different sweeping rates a = 1.0, 1.01, 1.05 and 
for each a we choose several magnetization quantities for 
comparison. From fig. 3, one can read that the spin tun- 
nelling happens mainly around B = regardless of the 



different quantities of magnetization, and we also see the 
tunnelling processes are very sensitive on the sweeping 
rates. Moreover we find the conserved magnetization m 
dramatically affects the tunnelling processes as well as 
the final tunnelling probability. For example, in the first 
row figures of fig. 3, with increasing the magnetization 
from to 0.3, we find that the occupation population of 
BEC in zero-spin component after a round sweeping of 
the external magnetic field changes from being enhanced 
to being quenched compared to its initial state. The in- 
fiuence of the magnetization parameter can be also seen 
from fig. 2, where the fluctuation on the tunnelling prob- 
ability is clearly suppressed by increasing the value of the 
magnetization. 

The crucial effect of the magnetization parameter on 
the spin tunnelling can be roughly understood from 
eqs.([Sl [7]). It shows that the variation of the population 
Sq is restricted by the conserved magnetization quantity, 
i.e., |1 - Sol > \m\- 

To explain why the spin tunnelling happens mainly 
around zero value of the magnetic field, i.e., i? = 0, we 
calculate the eigenvalues as well as the eigenstates of the 
system using the similar methods developed in our recent 
work [Ssjl . Solving the eigen equations of the system, we 
obtain the eigenvalues or eigenenergies. Fig. 4(a) plots 
them for linear case, i.e., c = 0. Fig. 4(b) plots them for 
nonlinear case, i.e., c 7^ 0. One can see, at c = 0, the 
system has three levels. They are ei = £'+ + c, £2 — 
Eq +c, £3 = +c, which cross around B = and corre- 
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FIG. 3: The change of so (Its initial is 0.5 at B ^ —00) with 
the sweeping magnetic field B for Rb at m — 0.0,0.1,0.3, 
a = 1.0, 1.01, 1.05 respectively, and c = -3.13 x 10"^* J. 
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FIG. 4: (Color on line) The eigenvalues of '^^ Rb atom con- 
densate, (a) is for c = and different m. (b) is for 
c = —3.13 X 10~^*J and differen m, where the triangle is 
£2a and the circle is £26- The inset shows £2a and £26 are non- 
degenerate, (c) is for c — —3.13 x 10~^* J and m = 0, where 
the dash line is Hmfi and the solid line is Hmf2- To clearly 
show the spin tunnelling around _B = in fig. 3, the unit of 
horizons of (a) and (b) is taken to be 10~^G, while it is G in 
(c). 



happens in a small regime around B = 0. By the way, 
from the above analysis we see that in this system the 
eigenstates do not correspond to the extreme values of 
energy, e.g., the state of sg = 0.5 in [13] is not an eigen- 
state but a state with an extreme energy. For the same m, 
tunnelling happens between these extreme energies, while 
the eigenstate on the eigenvalue (e.g., £26 in fig. 4(b)) is 
always not change. 
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spond to m = —1, 0, 1 respectively. When c ^ 0, the mid- 
level £2 is split into two levels, i.e., e2a = {E+ + E-) /2—c 
(the triangles in fig. 4(b)) and £2fc — Eq ~ c (the circles 
in fig. 4(b)) corresponding to m = {E+ — E-) /Ac and 
m = 0, respectively. e2a corresponds to the states of 
So = 0, while £26 corresponds to those of sq = 1. They 
are seemingly degenerate, while the inset graph in fig. 4(b) 
shows they are actually non-degenerate. Furthermore, 
we find an interesting phenomenon for this irregular sys- 
tem through investigating the extreme energies of the 
classical Hamiltonian. After adding the nonlinearity to 
the system, these extremes are different from the eigen- 
values for the same m. Taking m = as an example, 
in order to ensure the exact position of the extreme en- 
ergies, we use -ffm/ and obtain its two extreme values 
Hmfi = Eo, H^f2 = (5/2 + 6^ /16c + Eq + c. They are 
plotted in fig. 4(c). Due to these levels are very close 
around B = 0, tunnelling between these levels easily oc- 
curs. Moreover, our calculation reveals that these almost 
degenerate solutions are only emerging in the magnetic 
field of range [—0.16, 0.16] G. It means that the tun- 
nelling should mainly occur in this regime. The above 
analysis coincide with the jumping regime of sq in fig. 3. 
In this way, we explain why the spin tunnelling mainly 



FIG. 5: The relationship between the initial value of so at 
B ^ -oo, Sf^"'*^"', and its final value at B ^ +oo, sp'^"' , for 
a = 0.01, 0.5, 100, m = 0.0, 0.1, 0.3 and c = -3.13 x 10"^'' J. 

The above phenomena are general and is independent 
on the initial condition, namely Sq. To more clearly com- 
prehend this point of view, we investigate the relation- 
ship between the initial value of sq and s^"""'. In fig. 5, 
we take a = 0.01,0.5,100 as an example and calculate 
the relationship between initial sq and s^™"' for several 
TO. We see it is a smoothly diagonal line at a = 0.01, 100, 
which stands for no tunnelling. At a = 0.5, irregular tun- 
nelling occurs. For the same a, the larger m, the more 
smooth the line is, which indicates that to suppresses the 
irregular tunnelling. 

B. ^^Na Atom Condensate 

^■^iVa atom condensate has also been prepared in ex- 
periments, and its dynamics are also interesting. Dif- 
ferent to the interaction between ^^Na atoms is 
attractive. In this subsection, we will discus its spin tun- 
nelling. In our discussions, we first still take the initial 
probability of sq as 0.5 in order to compare with ^'^ Rb 
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atom condensate. Like the above subsection, our discus- 
sions are divided into three parts according to the values 
of the sweeping rate a. The main results are shown in 
fig. 6. We see the tendencies of tunnelling probability for 
"^^Na atoms are the same as those of ^^i?6 atoms. Two 
little differences are found by comparing the two atom 
condensates. 



in fig.8(b)(e)(h) as well as in fig.5(b)(e)(h). The differ- 
ence is that the amplitude of the fluctuation. This can be 
seen in fig. 6 and fig. 2. The suppression effect of the finite 
magnetization m on the irregularity of tunnelling prob- 
abilities for ^'^TVa system is less significant than that of 
system. When a is larger, i.e., in the sudden limit, 
fig.8(c)(f)(i) and fig.5(c)(f)(i) show identical smooth lines 
without any fluctuations. 
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FIG. 6: The tunnelling probability of Na for m = 
0.0,0.1,0.3 aX c — 3.13 x 10"'^'* J in the sweeping magnetic 
field. 

The first difference is the adiabatic range. For 237Va 
system, it is smaller than that of system. Because 
the frequencies of the instantaneous periodic orbits of 
'^'^Na atom system is smaller than those of conden- 
sate, the system needs a less sweeping rate to satisfy the 
adiabatic condition. This difference does not affect the 
fact of no tunnelling phenomenon at a ^ 0. Fig. 7 mag- 
nifies an adiabatically chaotic-like part of fig. 6(a) and 
shows the above phenomenon. When a — > oo and a is 
moderate, the tunnelling phenomena of '^^Na system are 
similar to the counterparts of atom system. 

The second difference is the effect of the conservation 
TO. Like atom condensate, the relationship between 
gF'inai ^Initial studied to morc clearly understand 
the irregular spin tunnelling of '^^Na atom condensate. 
Using the same values of a and to as that in atom 
system, fig. 8 shows the relationship. When the sweeping 
rate is small, comparing fig. 8 with fig. 5, we see the lines 
in fig. 5 are smoother than those in fig. 8 for a same a, such 
as fig. 5(d) and fig. 8(d). This indicates that the effect of 
the conservation to on the spin tunnelling of '^^Na sys- 
tem is less important. When a is moderate, finite value 
of TO suppresses the amplitude of tunnelling probability 



FIG. 7: The magnified part of the tunnelling probability of 
^"^Na at m = 0.0, c — 3.13 x lO"'^* J in the adiabatic sweeping 
magnetic field. 



IV. INTERPRETATION OF THE IRREGULAR 
SPIN TUNNELLING 

In this section, we achieve insight into the irregular 
spin tunnelling effect of spinor BECs with the phase 
space of the classical Hamiltonian He [13] ■ As is discussed 
above, the tunnelling phenomena of '^^Na and atom 
condensates have no essential difference. So we take 
system as an example to interpret the irregular spin tun- 
nelling observed in the above sections. 

Fig. 9 plots the phase space of Hamilitonian ([5]) for 
for different 6/c. In these phase space we can find two dif- 
ferent dynamical regions: (I) running phase region where 
the relative phase 9 varies monotonically in time; (II) os- 
cillation region where 9 oscillates in time around a fixed 
point. As S/c varies, the areas of these two regions change 
respectively. 

At first, S/c oo, all the trajectories are in running 
phase region. As B changes with a moderate sweeping 
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FIG. 8: The relationship between s^"'*^''' and sp'^"'' , for a = 
0.01, 0.5, 100, m = 0.0, 0.1, 0.3 and c = 3.13 x 10"^* J. 




FIG. 9; The phase space of Rb for different magnetic field 
B at c = -3.13 X lO'^^'j, m = 0.0. The unit of B is Gauss. 



rate a, we record sq, 9 after an interval of time and plot 
them in phase space (see in fig. 10). We see that main con- 
tribution to the spin tunnelling comes from the transition 
point where the trajectory passes from the running phase 
region to the oscillation region or vice versa. For differ- 
ent a, the transition point and the final equilibrium place 
are quite different, for example a = 1.0 and a = 1.01 in 
fig. 10. So the tunnelling probabilities are expected to 



be sensitive on the sweeping rates with showing irregular 
patterns observed in Fig. 5 and 8. 

For higher resolution of a, the seemingly chaotic tun- 
nelling probability is regular. In fig. 11 we plot the mag- 
nifying part of fig. 2 around a — 1.55. We find that 
around a precision a with high resolution the irregu- 
lar structure becomes the regular one having periodic 
structure. Set its period is ap. When the a is fixed, 
we find ap has following relation with the initial value 
of the magnetic field Bq: ap ^ l/|J5o|^. Fig. 12 shows 
the relationship between ap and Bq around a = 1.55 
and a = 2.55. They are 0.05/ |Bo|^ and 0.15/ \Bo\^ re- 
spectively. Furthermore, for a fixed Bq, the relation- 
ship between ap and a is ap ~ a^. Fig. 13 shows 
ap = 0.0004,0.0011,0.0022 around a = 1.55,2.55,3.55, 
respectively. We find 0.0004 : 0.0011 : 0.0022 is nicely 
equal to 1.55^ : 2.55^ : 3.55^. So, the period of regular 
structure near a sweeping rate is 

ap ^ oCk^ (8) 

\Bo\' 

From this formula, we see that a large Bq (used to cal- 
culate tunnelling probability in the above section) leads 
to a small ap. If the resolution of the sweeping rates is 
small compared to the above period, tunnelling probabil- 
ities are usually recorded in different period with random 
phase. So they will looks like chaotic. Only when the res- 
olution is high enough compared to the above period, the 
regular structure can be observed. 

Actually, for a moderate sweeping rate, the chaotic- 
like relationship between g^"-'*'"' and s^™°' in fig. 5 and 
8 is also due to the above reason, that is, when the res- 
olution of the initial values of s^"***'*' is increased, the 
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FIG. 11: The magnified part of the tunnelling probability 
of ^''Rb at m = 0.0, c = -3.13 x IQ-^^J in the moderate 
sweeping magnetic field. 



Sq ^ 1/ \Bq\. In fig. 15, this relation is confirmed by our 
numerical simulations even for different magnetization m 
and a same a. So we expect that the above inversely pro- 
portional relation is independent on m. 
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FIG. 12: The relationship between ap and Bo around a = 
1.55 and a = 2.55. The dots and the squares are numerical 
results. The solid line is ap = 0.05/ |Bo|^. The dash line is 
ap =0.15/|Bo|^ 



FIG. 13: The magnified part of the tunnelling probability 
of ^"^Rb at m = 0.0, c = -3.13 x 10"^* J around a = 
1.55, 2.55, 3.55 from up to down, respectively. 



V. CONCLUSION 

In conclusions, theoretically, we have investigated the 
tunnelling dynamics of a spin-1 Bose-Einstein condensate 
in a linearly sweeping magnetic field within a framework 
of mean field treatment. We focus on the two typical 
alkali Bose atoms ^"^Rb and '^^Na condensates and study 
their tunnelling dynamics according to different sweeping 
rates of external magnetic fields. We also investigate the 
effect of the conserved magnetization on the dynamics 
of the spin tunnelling. We hope our studies would bring 
more attention to spin tunnelling experimentally. 



observed irregular patterns observed are expected to dis- 
appear. Fig. 14 shows the magnifying part of fig. 5 around 
gi^nitiai ^ Q_4_ Fig.l4(b) and fig.l4(c) plot this kind of 
regular structure for different initial magnetic field. Set- 
ting the period of this regular structure as Sq, wc find 
it satisfy the following relationship between Sq and Bq: 
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